Space dependent Fermi velocity in strained graphene 
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We investigate some apparent discrepancies between two different models for curved graphene: 
the one based on tight binding and elasticity theory, and the covariant approach based on quantum 
field theory in curved space. We demonstrate that strained or corrugated samples will have a 
space dependent Fermi velocity in either approach that can affect the interpretation of local probes 
experiments in graphene. We also generalize the tight binding approach to general inhomogeneous 
strain and find a gauge field proportional to the derivative of the strain tensor that has the same 
form as the one obtained in the covariant approach. 
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As it is well known the low energy electronic excita- 
tions in graphene are modeled by a masslcss Dirac Hamil- 
tonian that is able to describe many spectroscopic and 
transport experiments with surprising accuracy The 
Fermi velocity is the only free parameter in the model 
and, as such, all the observable quantities depend criti- 
cally on it. It plays the same role as the effective mass 
in the usual Fermi liquid theory. Recent experiments J2 - 
Q have been able to measure the renormalization of the 
Fermi velocity predicted in [f| what means that the Fermi 
velocity is not a constant but increases as the energy is 
decreased near the Dirac point. In this work we will show 
that the corrugations and strain present in most of the 
samples give rise to a space-dependent Fermi velocity, a 
fact that can change the interpretation of some experi- 
ments. 

The presence of ripples in all the graphene samples and 
the influence of the lattice distortions on the electronic 
properties of the material are two of the most interesting 
aspects of graphene that still remain as open problems 
in the field. The most popular model proposed in the 
literature to study the issue is based on a combination 
of tight binding and elasticity theory [f| @] that induces 
"clastic" vector fields coupled to the electronic density 
(for an extensive review on the subject see Q and the 
references therein). This approach has given rise to the 
proposal of "strain engineering" [j| [l(J. An important 
highlight in the field has been the observation of a reor- 
ganization of the spectrum resembling Landau levels in 
strained graphene due to the effective magnetic fiel d fill 
whose existence was predicted theoretically in [To|, Il2|. 

An alternative model to explore the influence of corru- 
gations on the electronic properties of graphene is based 
on the formulation of quantum field theory (QFT) in 
curved space well known in cosmology and gravitation 
[l3l Il4j . This possibility, suggested for modeling spher- 
ical fullerenes in [la [lq | , was fully explored for general 



geometries m 




The QFT approach is rooted on 



the spinorial nature of the low energy electronic excita- 
tions of graphene and on the robustness of the Fermi 
points under deformations of the lattice [2l[ , and hence 



it should work well in the low energy scale. 

Since both models are natural and predictive, one 
should expect that they will provide the same results 
when applied to curved graphene samples with given 
shapes. Nevertheless there are a few discrepancies that 
are immediately apparent without entering into much de- 
tail. One is the prediction of the QFT approach of a 
space-dependent Fermi velocity put forward in [l9| . The 
origin of it lies on the vector indices of the Pauli matrices 
that made them space dependent when going to curved 
geometries. There is also a homogeneous contribution 
coming from volume effects. Another obvious discrep- 
ancy is related to in-plane distortions. These give rise in 
the elasticity approach to the coupling of the spinor to 
vector fields (named fictitious magnetic fields in the lit- 
erature) which in turn have observational consequences. 
Since in-plane distortions do not induce intrinsic curva- 
ture to the sample, the analogue of the fictitious vec- 
tor fields in the QFT approach, the spin connection, 
is zero. A third apparent discrepancy comes from the 
different symmetry of the vector fields associated to in- 
trinsically curved samples where the two approaches can 
be used. The best example is provided by the gaussian 
bump worked out in [l9| whose associated vector field is 
axially symmetric in the QFT scheme and has trigonal 
symmetry in the TB-elasticity models (22j . 

Although the tight binding derivation of the low energy 
effective Dirac Hamiltonian is the most popular proba- 
bly for historical reasons (23|, it has been known since 
1958 that the Dirac structure is more general and follows 
from the lattice symmetry and a low energy expansion 
p^ . New directions studying the effective low energy 
Hamiltonians for distorted lattices on symmetry grounds 
[25Tj29| paved the way to a better understanding of the 
correspondence between the lattice formulations and the 
QFT covariant approach. In the present work we par- 
tially follow this path and explore the correspondence 
between lattice and continuum formulations with special 
emphasis on the spacial dependence of the effective Fermi 
velocity. We will show that a space dependent Fermi ve- 
locity arises from the tight binding-elasticity approach 
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when going beyond the linear approximation and that 
the two formalisms can be compared if a metric coming 
from elasticity is chosen in the QFT approach. 

The QFT geometric description of curved 
graphene and the elasticity theory. 

The QFT geometric description of curved graphene has 
been discussed in detail in [30( ■ It is based on the stability 
of the Fermi points of the hexagonal lattice under mod- 
erate lattice distortions [2l| and on the subsequent de- 
scription of the low energy excitations around the Fermi 
points as masslcss Dirac fcrmions. Hence a natural way 
to incorporate the effect of the observed corrugations at 
low energies is to couple the Dirac equation to the given 
curved background. 

To make the connection between the two approaches 
explicit, we will work out the QFT Hamiltonian in a 
curved space arising from a metric related to the strain 
tensor by ay = rjij + 2uij [3lL [32| , where rjij is the fiat 
metric (the identity matrix), and the strain tensor is de- 



fined as 



djiii 



dihdjh), where Ui and h 



are in and out of plane displacements respectively. Since 
the metric g is already phenomenological, we could in- 
clude in its definition a material dependent parameter ft 
similar to the one obtained from the TB approach. We 
choose to leave its value to 1 for the time being. 

The Dirac Hamiltonian in a curved space described by 
a metric gij is given by 



H = i rf 2 x^fel(9, + l],)i 



(1) 



where ip = ^7°, 7°7 4 = <J l are the Pauli matrices, and a 
sum is implicit in repeated indices. The effect of the met- 
ric is encoded in the tetrads e* , the metric determinant 
yfg, and the spin connection f2,. Expanding these objects 
to first order in ity (see the Supplementary Material A) 
we obtain the Hamiltonian 



model. The general tight binding Hamiltonian is 

H = -Y] ts,n4 b s+L + cc -> ( 3 ) 



where x runs over the position of all unit cells and the 
three nearest neighbour vectors are defined as 

h=a{^,\) 8 2 = a{-^-,\) 53=0(0,-1), (4) 

with a the equilibrium nearest neighbour distance. 

In order to make an easier contact with previous works 
in the literature and to illustrate the method we will first 
consider the case of homogeneous strain tg n = t n and 
then generalize it to arbitrary inhomogencous strain. In 
the tight binding approximation, the only effect of strain 
is to modify the hopping integrals as the distance between 
atoms changes. Even if the strain does not depend on 
position, the three nearest neighbour hoppings t n may 
vary independently. To first order in the distance change 
Ait n , we may write 



t n = *o(l - /3Ait n ), 



(5) 



where (3 = |<91og£/<91oga| and to the equilibrium hop- 
ping. The relative distance change to first order in strain 
is 



2 "y 



(7- 



(6) 



Notice that when considering inhomogencous strain as 
done in the Supporting Information C this equation has 
to be completed with terms proportional to the derivative 
of the strain tensor. The lowest order in this case is 



6151 



5< S j S k 



2q2 -n,,,,,. 



(7) 



We now expand in Bloch waves 



H = if (Pxt/j* 1 (a a d a + Un<j a d a - a a u ab d b + (2) 
\a a {d a u ll - d l u ta ))ip. 

The first term in eq. ((2]) is the usual Dirac term. The 
next two terms are the space dependent Fermi velocity, 
and the last two are the corresponding geometric gauge 
field, as interpreted in ref. [lj|. Note that terms of the 
type i(uijdk + l/29feUy) are Hermitean by partial inte- 
gration: for every term in the Fermi velocity there must 
be a corresponding geometric gauge field that guarantees 
hcrmiticity. Also note the extra factor of i as compared 
to the minimal coupling of a U(l) gauge field, di + ieAi, 
which is Hermitean by itself. The effective Hamiltonian 
around the other Fermi point will have the same Fermi 
velocity but a minus sign in the gauge field couplings. 

Space dependent Fermi velocity from tight- 
binding. 

We now proceed to compute the effective Hamiltonian 
in the presence of strain directly from the tight binding 



kGBZ 



b s = 2^ e „ fi 



(8) 



kEBZ 



and consider momenta close to the Dirac points k = 
K + q, with K = (47r/(a3\/3), 0) (the result for the other 
Dirac point can obtained changing K — > —K through- 
out). The Hamiltonian in momentum space is 



H 







(9) 



Expanding to first order in q we obtain 



AK8 n 



ra=l 



iK6„ 







(l + i<r 3 q-6 n ). (10) 



We now use the following identity 




(11) 
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where a = (a x ,a y ) are the two Pauli matrices. Using 
eqs. (|5|) and ^ for t n the Hamiltonian is 



obtained in (|T6|) . The new term is a "geometric" gauge 
field given by 



H 



EM 1 



S„ud n J y-v^-Snj [l + i(T3q-5 r 

(12) 

We can now collect the different terms of this expression 
with the use of the identities given in the Supplementary 
Material B. Labeling the various terms by their order in 
the expansion in g, u we get 



H — Hg + H u + H q .u, 



with 



H q = VoaiQi 



2a 

VQ 



f3<T l K i jkekiUji, 



(13) 

(14) 
(15) 

(16) 



where we have defined vq = 3tna/2 and Kyk is the in- 
variant C3 tensor given in (|39| J8|. It is easy to see that, 
if needed, the expansion of the low energy hamiltonian 
in powers of qi and Uij can be done to any order. The 
first two contributions are well known: H q is the usual 
Dirac Hamiltonian, and H u is the standard strain in- 
duced gauge coupling Q: H u = v <JiAi with components 



A, 



A, 
2a 



"-yy) 



A, 



2a 



(-2u xy ). (17) 



The new contribution H u q is the main result of this work. 
This result is consistent with the one obtained from the 
symmetry analysis in [28l . |29| , and fixes the coefficients 
of the symmetry allowed terms in terms of a microscopic 
model. We will see its full significance in what follows. 

Inhomogeneous strain. Treating inhomogeneous 
strain is a delicate issue due to the lack of translation 
invariance of the system. The usual procedure [29[ con- 
sists in taking the homogeneous Hamiltonian in k space 
and go to real space by the replacement rule 



",j<lh -> i{u tj d k + l^dkUij), 



(18) 



which guarantees Hermiticity. Using this rule, the term 
([TBI) becomes 



■ Vo 



■13 



I- 

2 °% u jj 



what allows to write the total Hamiltonian as 
H = ir,,\r)n,i), + iv a t T t + V (TiAi, 



(19) 



(20) 



The field Ai is the one obtained in the standard approach 
(|17j). We also get the tensorial and space dependent 
Fermi velocity 



Vij 



^(2 Uij 



VijUkk) 



(21) 



r ' = 4 ( djUij + \ d%u i 



(22) 



which was obtained in the covariant approach ([2j . Being 
proportional to the derivative of the strain tensor it only 
appears in the case of inhomogeneous strain. All these 
terms are also found in the symmetry analysis 0, [29[ . 

As discussed in [33|, the former procedure can miss 
some terms. We have checked that this is the complete 
result to this order in derivatives by directly performing 
the Fourier transform of the real space Hamiltonian in 
the first stage of the tight binding for the general case of 
inhomogeneous strain. The details can be found in the 
Supporting Information C. 

Discussion and future. 

The purpose of this work was to discuss the equiva- 
lence of the tight binding and the QFT approaches in the 
description of curved or strained graphene with special 
emphasis on the space dependent Fermi velocity. The 
main result is that the two approaches give rise to the 
same type of terms for the spacial dependent Fermi ve- 
locity although with different numerical values of the co- 
efficients. The symmetry approach gives rise to the same 
type of terms with independent - and undetermined - co- 
efficients whose value has to be fixed by the model. We 
have also derived the " geometric" gauge field in the tight 
binding approach for the case of inhomogeneous strain. 

Referring to the QFT approach, we have seen that the 
definition of the metric of the curved space in terms of 
the strain tensor allows to get physical effects from the 
covariant formalism even in the case of having only in- 
plane distortions. 

An interesting related issue concerns whether new 
terms can appear in the effective tight binding Hamil- 
tonian independently of the change of the hoppings, pro- 
duced just by changes in the relative positions of the 
atoms. If one assumes (3 = 0, the discrete TB Hamilto- 
nian in that case is just 



(23) 



with i,j running through nearest neighbour atoms. It is 
easy to see that moving the atoms form their positions 
does not change the Hamiltonian at all. i,j are just la- 
bels numbering the atoms, and need not refer to physical 
position in any sense. Hence the energies and eigenfunc- 
tions of the system are not modified by this strain, even if 
inhomogeneous (for /3=0). However, this does not imply 
that there will not be observable consequences. The fact 
that the atoms are in different positions does change the 
way in which external position-dependent probes see the 
system. In short, when the atoms are displaced, the dis- 
crete label i maps to a physical position in the strained 
frame. If wc describe the physics in the lab frame so 
as to be able to couple an external field, vectors in the 
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strained frame have to be rotated to the lab frame. This 
gives rise to a /3 independent contribution of the type 
obtained recently in [34| and will be discussed in detail 
elsewhere. 

Although the space dependent Fermi velocity is ob- 
tained at higher order in a tight binding expansion, its 
presence has important physical consequences and it can 
not be obviated. The important issue in this discussion 
is the influence of the various factors entering the effec- 
tive Hamiltonian on the observable physical quantities. 
In particular in fl9j it was shown that in the QFT ap- 
proach the local density of states (LDOS) is not affected 
by the vector fields to first order in perturbation theory. 
The oscillations in the LDOS obtained in that work came 
from the combination of variable Fermi velocity and vol- 
ume effects. Similarly it has been shown in [35j that long 
range correlated vector fields do not alter the minimal 
conductivity of graphene that is in turn severely changed 
by long range correlated disorder in the form of a ran- 
dom distribution of Fermi velocity (3(| . We note also that 
while this paper was completed we learned on an exper- 



iment that points to the observation of a 5-10% spatial 
fluctuation of the Fermi velocity in samples on SiC>2 ■ 

As we have seen, the standard QFT approach is rooted 
in considering the spinor describing the low energy elec- 
tronic excitations of graphene as a covariant spinor under 
a geometric point of view. Since the time coordinate re- 
mains "flat" , the Lorentz symmetry is reduced to trans- 
lation and rotations. This approach has been pushed 
forward in (37j considering graphene as a QFT in curved 
spacetime what gives rise to very interesting consequences 
as the possibility to observe a Hawking-Unruh tempera- 
ture in the curved graphene samples. 
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I. SUPPORTING INFORMATION 

A. The QFT geometric description of curved 
graphene and the elasticity theory. — 

In this section we will work out the various geometric 
objects entering in the curved space Hamiltonian in eq. 
(1) of the main text. 

The volume element will be given by Jg = detgij = 
1 + uu. The next geometric object needed is the tetrad, 
which satisfies 



— ^j^ab — &i 7 



The easiest choice is 

e-ai = 8 ai + S J a Ui 

el 



'ij • 

J a — u a u aj u > 

Next we need the affinc connection 



(24) 

(25) 
(26) 



r fe _ 1 n ikr d 9ji , d 9u_ dgjj 1 . , 

which to first order is 

T% = 5 lk {d iUjl + djuu - d lUij }, (28) 
Next we need the spin connection, which is obtained from 



1 



1 



= llalbe a 3 g^ \e\ = --a a a b e^g^\7 ie l (29) 



with 



To first order the spin connection is 



1 



= -^o a o h {d a u ih - d b u ia ). 



(30) 



(31) 



With aid of the identity 

b c a cbc , c cba b rac , ■ abc /on\ 

a a a = a o + a o — a o + te , (32) 
we can compute 

a a e l a a l = -]^(Ta{d l u la - d a u u ). (33) 
The Dirac Hamiltonian to first order in u finally reads 

H = if d 2 x^{a a d a + Uua a d a - o a u a bdb + (34) 
\(j a {d a un - diU ia ))ip. 

B. Some tensor identities 

The following identities are used in the tight binding 
derivation: 



^E4# = ^'> ( 36 ) 

n=l 

i^#J = -^, (37) 
n=l 

^ E Wn 5 n = l(v ij V H + rf k rf l + V il rf k ). (38) 

71=1 

The object K l ^ k in the second equation in ([55} is an 
invariant tensor under the discrete C3 rotations of the 
lattice given by 

3 

R abc = J- d a d b d c = d a d b d c + jajbjc + jajb^ (39) 

n=l 

with di the nearest neighbour vectors. It is easy to see 
that the only non-zero components of it are 

111 1^122 7^212 K--221 



K 111 = -K lAi = —K 



K^ 1 = 1. (40) 



C. Inhomogeneous strain with explicit Fourier 
transform. 



We start with the standard tight binding Hamiltonian 

# = - E **»4w„ + hc -> ( 41 ) 

.r. n 

where x runs over the position of all unit cells and the 
three nearest neighbour vectors are defined as 

£=a(^,l) 5 2 = a(-^-,±) S 3 = a{0,-1), (42) 

with a the equilibrium nearest neighbour distance. We 
now replace a x and b x by their fourier expansions and 
get 

# = -£*,»£*,***,» (43) 

Using the expressions Uk = -k J2 X e~ %kx u x , and 
J2 X e lkx = NSk,o, H becomes 

# = -£*,»£*,*'**.» ( 44 ) 

[e-' li3 ale lk ' {s+ ^b k , + e^*' e ife a fe ) . 

To write this as a matrix equation we have to relabel 
k <H- k' in the first term, so that we get 



h = -n J2 (4»*b)t*-k,n 



fe.fe' ,11 



MTi<T3 = £ij<Jj 



(35) 



-iWS n J V bk 

(45) 
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Now wc proceed as usual, expanding k = K + q and 
P = K + (f. We get 

(46) 



To expand in q, q' we rewrite this as 



H = -N V fat, fc n t , / , + f , 



and using the identity 

/ q e -iK-S„ \ g.S 



iK-S„ 







i -<x 3 , (48) 



H = — ^ Cq>tq'-q,n(i<T3v8 n ) (l - ZtT 3 ( g "^ g )S n + )8 n \ C q , (49) 



we get 

where Ck = (a^, 6/.). The modified hopping can now be expressed as the Fourier transform of eqs. (5), (7) of the main 
text as: 

/ s n s n <p s j s k \ 

tq'-q,n = h ( fy.g " P U q'-q,ij^T ~ 1 ~ <l)kUq>-q,ij " j . (50) 

It is easy to see that the free Dirac term (linear in q and zeroth order in u) and the standard TB gauge field A (linear 
in u and k=K) are recovered from this expression. The new gauge field comes from the crossed term linear in u and 
<?: 

H u , q = 3h ^ N ]T c\, u q >- qiij i(X3a k (-ia 3 (qi + q[) + i(© - q'i))(v ij V M + rf k rf l + rf l rf k )c q 
q,q' 

+ E C S' ff 3 - q[) (v li V kj + '/V + rf k ^)c q . (51) 



r 



The terms with q— q' are equal and with opposite signs, with 
so they cancel each other. The fourier transform back to 
real space is done with the prescription 



^ + (52) /3 / 1 . 

T, = 7 djUij + . 54 



obtained from the definition of the Fourier transform. 
This allows to write the crossed term in the form 

H u ,q = i-^-( 2u i3 + VijUkk)°idj + iv <JiTi, (53) 



